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Abstract: We couple fermion fields in the adjoint representation (gluinos) to the SU(2) 
gauge field of unit charge calorons defined on x Si. We compute corresponding zero- 
modes of the Dirac equation. These are relevant in semiclassical studies of A/" = 1 Super- 
symmetric Yang-Mills theory. Our formulas, show that, up to a term proportional to the 
vector potential, the modes can be constructed by different linear combinations of two 
contributions adding up to the total caloron field strength. 
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1. Introduction 

In recent years a deep link between monopoles and instantons has been established at 
finite temperatm'e It was for long known that the finite temperature instanton, 

the Harrington-Shepard caloron becomes for large scale parameter a BPS monopole 
0. This happens when the value of the constant Polyakov loop at infinity, the holon- 
omy, approaches unity. The surprise came with the discovery that non-trivial holonomy 
calorons when squeezed in the time direction reveal to be composed of N, for SU (N), BPS 
monopoles. — 1 of these constituent monopoles are massless for the HS caloron, but 
in general all of them get non-zero masses with values related to the eigenvalues of the 
Polyakov loop at infinity. The idea of the composite nature of instantons, with instanton 
quarks or merons Q as constituents, has been on the basis of several semiclassical pro- 
posals to address the confinement problem in QCD. Isolated fractional instantons (twisted 



instantons) I IC] can be obtained by using tori with twisted boundary conditions []ll|. By 
replicating the tori one can obtain classical configurations in a periodic box, where the 
action density is clustered into lumps of of topological charge. These structures were 
observed in lattice generated ensembles at zero temperature and were argued to be rele- 
vant for QCD confinement |12]. Non-trivial holonomy calorons also exhibit explicitly this 



composite nature as far as the separation between constituents stays larger that their size. 
Otherwise they merge in an undissociated instanton. Triggered by this result, quite a num- 
ber of more recent lattice analysis have identified the presence of constituent monopoles at 
temperatures below but close to the deconfinement phase transition in Monte-Carlo gen- 
erated configurations ||l^, 0]. The possible relevance of the instanton-monopole link for 



- 1 - 



QCD dynamics is an open issue, nevertheless constituent monopoles have already shown 
their usefulness in a different context, in particular for calculations of the gluino condensate 



in 4D N = \ supersymmetric Yang-Mills theory |15]. The contribution of the constituent 
monopoles seems essential there to bring to agreement strong and weak-coupling instanton 
calculations of the condensate. 

The existence of an analytic expression for the non-trivial holonomy calorons allowed a 
subsequent analytic calculation of the zero-modes of the Dirac equation, in the fundamental 



representation, in this background fieldlTq, 17]. For large constituent separation the modes 



are entirely supported on just one of the monopoles, jumping from one to other as we 
change the periodicity condition in the time-like direction imposed on the solution. This 
knowledge has proven useful in interpreting the results of several lattice studies which 
employ low-lying eigenstates of the Dirac operator to trace topological structures present 



in gauge field configuration ensembles |13], |18| 



The present paper is devoted to the derivation of the analytic expression and proper- 
ties of the zero-modes of the Dirac equation in the adjoint representation for Q=l SU(2) 
calorons. These are relevant objects in the study of the semiclassical behaviour of 4D A/" = 1 
supersymmetric Yang-Mills theory compactified in R"^ x S^. They are directly related to 
estimates of the gluino-condensate [21|. By now this has been studied |15| only in the 



limit of large constituent separation, in which the SU(2) caloron degenerates into two BPS 
monopoles. Our expressions for the four adjoint modes, valid for any separation, indeed 
tend in this limit to two pairs, one attached to each constituent monopole reproducing the 
well known adjoint zero modes of BPS monopoles 0]. 

Our work is also useful within the previously mentioned spirit of using modes of the 
Dirac equation as probes of gauge field structure, an approach that has become very popu- 



lar (see for instance ||2^- [^) since the discovery of lattice Dirac operators |2g] that possess 



exact index theorems. The usefulness of adjoint modes in this respect has been recently ad- 



vocated ||2^. The idea is based upon the so called supersymmetric modes, having densities 
that match the action density profile but are less sensitive to ultraviolet fluctuations. 

The paper is organised as follows. In section 2 we describe the general properties of 
adjoint zero-modes, such as their occurrence in pairs, related by Euclidean CP transfor- 
mations. Then, we give the main formula for the modes within the ADHM formulation, 
and apply it to the Q = 1 SU(2) caloron case. This provides two pairs of zero-modes. One 
pair is given by the supersymmetric zero modes, which are proportional to the gauge field 
strength itself. The remaining pair of adjoint zero modes is studied. In addition to the 
analytical expression (details of its derivation are given in the Appendix), we display its 
density profile in some representative cases. In section 3 we show how the solutions behave 
in certain limits and how they interpolate between modes of BPS monopoles and those of 
instantons. We end up with a summary of the results and a list of possible extensions and 
applications. 

2. Formalism 

The adjoint zero modes ^%{x) are solutions of the Euclidean 4-dimensional massless co- 
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variant Dirac equation in the adjoint representation of the gauge group: 



^^ = (2.1) 

In this paper we will analyse the simplest case given by group SU(2). Then the colour 
index a takes the values 1,2,3, while the spinorial index a takes four values. Action by 
75 maps zero-modes into other ones. It is convenient then to combine zero modes into 
eigenstates of 75 with eigenvalue ±1, known as left and right-handed modes respectively. 
In Weyl's representation of the Dirac matrices, the left and right handed modes reduce to 
two-component spinors ip± satisfying 

= (2.2) 
Dip^ = (2.3) 

where D = cr^D^ and D = Wf^D/^i. The Weyl matrices are given by a^j = (I, —ir), in terms 
of the Pauli matrices Tj (the matrix is the hermitian conjugate of cr^). If the gauge field 
is self-dual, then Eq. |2]^ implies 

D^ij_ = (2.4) 

for all This can easily be shown to imply that the gauge-invariant density ['(/'-(a;)!^ must 
be constant (x- independent). For non-compact space-times these are non-normalizable 
solutions. 



Focusing now on left-handed zero modes (solutions of Eq. 2.3) we point out that the 
space of solutions is always even-dimensional. This follows from euclidean CP invariance 
mapping one solution into other 

^/J+^r+ = -ir2^l (2.5) 

In the previous formula -0+ stands for the complex conjugate spinor and the matrix T2 acts 
on the 2-spinor indices. Furthermore, for self-dual gauge fields one can establish a one to 
one correspondence between self-dual deformations (5A^ of the gauge field and left-handed 
zero modes. Given a deformation, one must first transform it to the background Lorentz 
gauge 

D^,5A^, = (2.6) 

and then one can show that 

= 6A^a^V (2.7) 

is a zero mode for any constant 2-spinor V . 

If there are isommetries of the problem which do not leave the solution invariant, 
the corresponding deformations are associated to specific zero-modes. In particular, the 
super- symmetric zero-modes arising for 5A^ = F^q, are associated with translation symme- 
try. In general, the space of self-dual connections is continuous and can be parameterised 
in terms of a set of real parameters (moduli). Variations with respect to these moduli 
parameters(tangent vectors) give rise to adjoint zero-modes. 
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Self-dual gauge fields with topological charge Q on the sphere 5*4 or on i?^ (with 
finite action) can be constructed by an algebraic procedure known as the ADHM con- 
struction [^8|. The fields are written in terms of g, a Q dimensional column vector of 
quaternions, and A, a. Q x Q matrix of quaternions satisfying certain conditions, to be 
specified later. In what follows, we will identify quaternions with the space of two by two 
matrices which are real linear combinations of the Weyl matrices a^. In particular, one 
can form the quaternion x = Xfj_a^ and its adjoint x = Xf^Wfj,. The self-duality condition 
amounts to the requirement that the matrix R: 

R = q0q^ + (A^ - x){A- x) (2.8) 

is real and invertible. 

The self-dual deformations are then associated with variations of the ADHM data q and 
A [^9|. Using the previously mentioned relation between adjoint modes and deformations, 
one obtains the formula for the modes in terms of variations of the parameters 5q and 5A: 

SAf, = ^{6q^ - u^6A)af'duj + h.c. (2.9) 

where we have introduced the x-dependent vectors of quaternions u and uj defined by the 
relations 

u = R-^q (2.10) 

and 

u = F{A-x)uj (2.11) 
where F = 1 + li^n is a real function. The symbol d in Eq. 2.£ stands for the contraction 



The condition that the variation 5Af^ given in Eq. 2.S satisfies the required covariant 
background gauge condition is 

Re (a^ 5 a -5A^ A + q{6q)^ - {5q)q^) = (2.12) 



where Re stands for the real part of the quaternion. This condition will be seen to hold 
for our formulas. Its interpretation will also become more clear later. 

Now we will proceed to particularise to the case of the Q=l caloron. This is a self- 
dual configuration in R^ x Si. At infinity the time- like Polyakov loop (the holonomy) is 
non-trivial. This is determined by the parameter 6i such that the trace of the Polyakov 
loop at spatial infinity tends to cos(27r(5i). 

For a given holonomy and a fixed period in time /3 (which we will henceforth fix to 
1), solutions depend on the following parameters: the position of the center of mass of 
the caloron Xqm, the size parameter p and an SU(2) colour orientation. For medium and 
large values of p (compared with the time-period /? = 1) the action density of the solution 
appears as a superposition of two lumps, named constituent monopoles in Ref. Q. The 
distance between the lumps approaches irp^ and the total masses Ma = ^Trnia/ are given 
by the holonomy as follows: 

nil = 47r5i ; m2 = 47r(52 = 2tt — mi (2-13) 
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The shapes of these monopoles tend, as the distance is increased, to that of BPS monopoles, 
having a non-abehan core which is exponentially localised and an abelian powerlike fall-off 
at large distances. 

One can construct the caloron solution by an infinite dimensional generalisation of the 
ADHM construction (strictly speaking a Nahm ||3^ transform). This is the procedure 
that we will follow here, allowing us to extend the ADHM formulas for the adjoint modes 
to this case. If we regard the caloron solution as a solution in R'^ which is periodic in 
time, the topological charge would now become infinite. Thus, q and A become an infinite 
dimensional vector and matrix respectively. The discrete index can be interpreted as the 
Fourier mode of a periodic function of one variable z (with period 1). Thus q{z) is a 
distribution and A a linear operator in this space. One can use translations, rotations and 
gauge transformations to bring q{z) and A to the form ||l|: 

q{z) = p{5{z - 6i)P+ + 6{z + Ji)P_) (2.14) 
A{z) = - iX'rxiiz) - iX^TX2{z) (2.15) 

where P± = 5{z) are periodic delta functions and Xa{z) are characteristic functions of 
the intervals la (taking the value 1 in the interval and zero elsewhere). The intervals Ii = 
[—5i , 5i\ and I2 = , 1 — denote complementary regions of length nia/ {mi + 771,2) within 
one period in z. Finally, the vectors X° can be interpreted as denoting the spatial locations 
of the constituent monopoles. We have used the translation and rotation symmetry to place 
them along the z axis and to locate their center of mass at the origin (miX^ +m2X'^ = 0). 
In addition, their separation is fixed by p: 

=V(0,0,1) (2.16) 

This information allows to determine X"" uniquely. As mentioned all Q = 1 caloron so- 
lutions can be obtained from these formulas by applying euclidean and gauge transfor- 
mations. Furthermore, one can easily restore of arbitrary time period by multiplying all 
length parameters by (5 (and masses by 1//3) . 

Within the Nahm transform philosophy, the quantity A can be identified with the 
covariant derivative (divided by 27ri) of an abelian gauge potential A^ over a 4-d torus 
which has been shrunk to a circle, whose coordinate is labelled by z. The remaining 
(spatial) coordinates have dropped as arguments, but the vector potential field still keeps 
the vector index. In our case only the third component is non-zero, and is given by 

I3 = -2tt{xIxi{z) + Xlx2{z)) (2.17) 

This implies that the corresponding magnetic field vanishes and the electric field is a delta 
function over z = ib(5i. 

In conclusion, to obtain the expression for the adjoint zero modes for our case, one 



has only to substitute the expression of the variations 5q and 5 A in the formula Eq. 2^ 
The variations are associated to the parameters of which the caloron field depends. On 
one hand we have the coordinates of the center of mass of the constituent monopoles. This 
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will give rise to the supersymmetric modes, which are always associated to translational 
symmetry. The corresponding variations can be obtained straightforwardly, giving 
and 6^^^ A = I. Substituting into the formula one obtains 6^^^Ao = and 

<5(°Ui = Ei (2.18) 

where Ei are the components of the electric (or magnetic) field strengths of the caloron. As 
anticipated this is the expression of the supersymmetric mode, having density proportional 
to the action density of the caloron field. 

The index theorem suggests that we should find 4 independent solutions {2NQ). As 
mentioned previously they come in CP-pairs. Each pair is associated to 4 real variations 
of the Nahm data. This is exemplified with the supersymmetric modes, for which there is 
a single pair associated to the 4d center of mass variations. Therefore, one must still find a 
new independent CP-pair. As we will see, one can obtain one such mode by varying with 
respect to the parameter p. Acting with the operator on our expressions of q and A, 
we obtain 6^^\ = q and 

5(1)1 = -2iXlT3Xi{z) - 2iXlT3X2{z) = ip\^{m2Xi{z) - miX2{z)) (2.19) 



For these variations q{5'^'^'> q)"^ — {5^^^q)q'^ vanishes, and Eq. 2.12| amounts to the requirement 



that the Nahm-dual field Afj_{z) satisfies the covariant background gauge condition too. 
This follows trivially, since (5(iUo = 0. In summary, substitution of these variations into 
Eq. |2.9| provides a self-dual deformation in the covariant background Lorentz gauge. The 
same applies to the supersymmetric zero-mode for which 5^^^q = and S^'^'Aq is constant. 

To give a simple expression of the result it is convenient to separate the contribution 
of 5q^ present only in the non-super symmetric zero-mode case from the other one. Further- 
more, we point out that the former becomes proportional to the caloron vector potential 
itself: . 2 

^-q^a^^duj + h. c. = jA^ (2.20) 

To give an expression of the second term one must realize that both u and uo can be 
written as the sum of two contributions, u = Ua and oj = Yla such that each piece 
is proportional to the characteristic function Xa{z) of each of the intervals. Finally, we can 
collect the formula for both sets of modes as follows 

5^^^A, = El + El (2.21) 
5^^^ A, = -jA^ + pWa^{^2Ei - miEl) (2.22) 

where rja^ is 't Hooft symbol and we have defined 

K = 'i^uiaadiVa + h. c. (2.23) 



which by virtue of Eq. 2.21 can be regarded as the contribution of constituent monopole 
a to the caloron field strength. This formula is quite appealing since it suggests that the 
two modes are simply given by different linear combinations of the field produced by each 
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constituent monopole. This is modified by the presence of the rjj^ in the expression for 
the non-supersymmetric mode. Notice, however, that each mode belongs to a complex two 
dimensional space of modes generated by the two elements of a CP-pair. In particular, this 
means that one can transform the gauge variations as follows 

= rj''/5A^ (2.24) 

for any value of a. Using this and allowing for a different normalisation one can recast the 
formula for the non-supersymmetric mode as follows: 

6'('^A^ = -nl^-^A^ + {m^El - mii?2)/2vr (2.25) 

In the appendix we compute the expressions for the functions Ua, ooa and with them 
we compute E%. The final expression Eq. ^.39| is a sum of two terms. The first one is 
given by the field produced by a BPS monopole, of mass Ma, located at the position of 
the corresponding constituent monopole, gauge rotated and weighted by an x-dependent 
scalar function \a{x). In the next section we will analyse how these functions behave in 
different limits and reduce to the formulas for monopoles and instantons. 

All of the expressions are dependent on a 2 x 2 matrix which is additive with respect 
to the contributions of the intervening monopoles. Its inverse, labelled V , appears in the 
formulas and provides the main effect of one constituent monopole over the other. 




Figure 1: Supersymmetric (left) and non-supersymmetric (right) zero-mode densities for 5i — 0.2 
and p = 0.1. Monopoles are localised on the z axis at x = y = t = 0. Densities are plotted in the 
X — z plane keeping t = and y = 0. The lengths of the x and z axes are scaled by 4.6p. 

As illustration of the properties of the zero modes we show in Figs. (|l|)-(0) the densities 
of both supersymmetric and non supersymmetric zero modes for Si = 0.2 and several values 
of the scale parameter p. Densities are plotted in the x — z plane keeping t = and y = 0. 
For small p the supersymmetric zero-mode reproduces the characteristic single-instanton 
shape. As p increases the caloron dissociates into two constituent monopoles which tend. 
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at large p, to two BPS monopoles. The non-supersymmetric mode has, at smah p, a 
symmetric ring structure also characteristic of a normal instanton, going through zero at 
the center of mass of the caloron. The ring gets distorted as p increases and dissociates 
for even larger p into the two constituent BPS monopoles. This behaviour matches the 
one obtained analytically in the p ^ and p ^ oo limits, which are described in the next 
section. 




Figure 2: The same as in Fig. pi but for p ~ 0.3. 




2.04 2.04 



Figure 3: The same as in Fig. y but for p = 0.8. 



3. Properties of the solutions 

In this section, we will clarify the spatial structure of the modes by studying their be- 
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haviour in certain limits. Since calorons interpolate between monopoles, instantons and 
HS calorons as one moves along the moduli space, a similar phenomenon is expected for 
the corresponding gluino zero- modes. 



3.1 p — > CO limit 

There are several length scales involved in this problem {m~^, the time period and the 
distance among constituents vrp^), so it is important to clarify precisely in what regime we 
are thinking of. We will consider the simplest situation in which irp'^ is much larger than 
all other length scales in the problem. In addition, we will focus on the behaviour close to 
one of the monopoles ri « np'^. Then, for the other monopole one has r2 ~ vrp^. 

All our expressions are built in terms of V and Ua (see Eqs. 5.28 and f).29 in the 



appendix), so we will first analyse the behaviour of these quantities as a function of p. 
The quantity U2 behaves for large distances as r2/27r — > /0^/2. Thus, V becomes order 
p~'^. In order to compute the leading behaviour of F in this limit one has to keep the first 
correction as well. This is obtained by expanding 

r2 = ||ri - 7rp2fc|| = np^ - (xs - X^) + . . . (3.1) 

In this way we arrive at 

V = ^il--,{U,-^-^^^^)) (3.2) 

and from here we see that F is order p'^. From these results we conclude that W (Eq. [5?2^ ) 
is order p~^, and dW of order /)~^ while the derivatives of Ua are order 1. Then it is easy to 
see that, of the different terms contributing to the zero-modes, only the term proportional 
to the a = 1 BPS monopole field is leading order. The latter is easily computable by 



realizing that the quantity inside parenthesis in Eq. 5.37 now becomes precisely 2/{Fp^). 
This leads to Ai = 1 and A2 = 0. The conclusion is that for large separations and close to 
the center of the ath constituent monopole, the quantity is, after a performing a gauge 
rotation, simply given by the electric field of a single BPS monopole: 

E^^^{x-X''-ma) (3.3) 

Now we recall that the integral of the BPS monopole electric field square over space equals 
ATim. We reproduce, in this limit, the well-known fact that the contribution of each con- 
stituent monopole to the energy will then become, at sufficiently large separation, propor- 
tional to mi and m2 respectively. Adding the two contributions to get the supersymmetric 
mode gives a total integral of Svr^ as expected for the caloron. 

Considering now the non-supersymmetric mode and taking into account that in Eq. 2.22| 



the first term is sub leading in this limit, we conclude that the contribution of the ath con- 
stituent monopole to the total energy is given by 47r p^mfml/ma. This is consistent with 
the known result about the integrated densities which appear when computing the 
metric of the caloron moduli space. 
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3.2 p — >0 limit 

To be precise the limit is defined by the requirement that maVa << 1 and the distances 
are of order p. This implies that in the first approximation one can neglect the separation 
between the constituent monopoles (vrp^ << r-j). As usual we will start by computing Ua 
in this limit: , 



eaXo 



Ua={ ^T" ) (3.4) 

where = Xq + is the square of the 4d euclidean distance to the origin. Notice that 
the 11 element of this matrix is negligible with respect to the 22 element, but has to be 
kept to be able to compute the inverse matrix. From here we conclude that to leading 
order V = 2Pii/(x^ + p^), where Pu is the matrix (projector) whose only non- vanishing 
component is 11, equal to 1. This projector appears in all our formulas and simplifies 
all vectors and matrices into single component. Using this it is easy to estimate the 
contribution of both terms in the expression for Ef^, Eq. |]3|. The second term turns out 
to be subleading with respect to the first one. Using d{Ui + U2) = x the second term in 
this equation gives: 

^-^^^ + /..c. (3.5) 

27r (x2 + p2)2 3,2 

which is proportional (ma/(27r)) to the gauge field of an instanton in a certain gauge. 
Hence, adding the a = 1 and a = 2 contributions, we conclude that the super symmetric 
zero mode reduces in this limit to the one of an BRST instanton. Notice however that the 



combination of Eq. |3.5| entering in the non-super symmetric zero mode vanishes. 

For the case of the non-supersymmetric zero mode the leading term turns out to be 
the one proportional to the vector potential. For it we obtain 

F 2[x^ + 

Notice that this is precisely proportional to the well-known non-supersymmetric adjoint 
mode for the BRST instanton. 

In summary, in this limit both adjoint modes tend to the corresponding ones for a 
BRST instanton. 



4. Conclusions 

In this paper we have derived, following the ADHM formalism, the analytic expressions 
for the gluino zero-modes of the Dirac operator in the background of topological charge 1 
calorons with non-trivial holonomy. These gluino zero-modes are relevant for semiclassical 
studies of 4D Super-symmetric Yang-Mills theories. They can also turn out to be useful 
for analysing the structure and topological content of the QCD vacuum. In particular they 
can give a handle in the identification of constituent monopoles inside instantons, a subject 
that has recently received much attention [|T^] -|pO[|. 

For Q = 1 there are four linearly independent zero modes which come in pairs related 
by Euclidean CP transformations. Two of them correspond to the super-symmetric zero 
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modes and share the property that their density exactly reproduces the action density of the 
caloron. The other two modes show a quite distinct behaviour. For small scale parameter 
p they have the same ring structure as for trivial calorons or instantons as exhibited in Fig. 
|l]. As we increase p the non-supersymmetric zero- mode rings get distorted (Fig. ^). For 
much larger p both zero mode densities dissociate into two structures having the density 
profiles of BPS monopoles of masses Mi and M2 (Fig. For the supersymmetric case 
the integrated contribution of each monopole to the energy is proportional to its mass. 
However, for the non supersymmetric modes these contributions are inversely proportional 
to their masses. By taking appropriate linear combinations it is possible to construct zero 
modes which, for intermediate to large scale parameters, single out only one of the two 
constituent monopoles. 

All these features can certainly help in the identification of constituent monopoles 
in lattice generated ensembles. This identification has been up to now performed solely 
on the basis of the action density and the properties of fundamental zero-modes |T3| , p!8| , 
p^ , becoming particularly complicated for low temperatures [20|. Adjoint zero modes 
provide an additional tool for this analysis. The supersymmetric zero mode density gives 
an estimate of the action density of the gauge field itself which is less sensitive to ultraviolet 
divergences. 

Our present results can be extended to the case of SU(N) Q=l calorons and to higher 
charge calorons. Furthermore, it is also possible to use our techniques to construct gluino 
zero-modes which are anti-periodic in time, directly relevant for finite temperature M = 1 
supersymmetry [32]. 
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5. Appendix 

In this appendix we will present the derivation the analytic formulas of the zero modes 
given in the text. All interesting quantities are expressed in terms of the functions u and u; 
(see Eqs. 2.11- 2.10| ), and their integrals over z. These functions belong to a 4 dimensional 



vector space over the field of quaternions. The space is the direct sum of the spaces of 
functions annihilated (up to delta functions) by the operator M = A — x and its adjoint. 
To obtain a basis of this space one must consider the functions 

^>^+\z,x,5)=x{-5,5)(^^-^' (5.1) 
*(-)(z,x,<5)=x(-<5,<5)e^'-^^ (5.2) 
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They are to be taken as periodic over z ( and xq) with unit period. The function x{~^i ^) 
is the characteristic function over the interval. In terms of these functions we can compute 
the 4 elements of the basis of the space: 

= (5.3) 

47r 

The index a takes two values (1 and 2), which can be thought as labelling the two con- 
stituent monopoles. Indeed, X° stands for the location and ttt-q = ^nba for the mass (for 
= 47r) of each of the monopoles. The remaining coefficients are = and = ^. 
To derive the formulas one needs to know how the operators M and act on them. 
These are given by 

M^^t = ^ ((e? + ie2n^ - 6i) + (e? - ie^^)6{z + <5i)) (5.4) 
M^a = ^ m + ie^2nz - 6i) + (e? - ie^,)d{z + 6i)) (5.5) 
(M - M^)^^ = 2i{x - X")r (5.6) 

The quantities e" and 62 are simple quaternionic functions of x defined by the value of the 
basis functions at the extreme of the interval: 

lim ^+(z) = ieaiet + ie|) = {ef + ie'^^) (5.7) 

where ei = — 1 and e2 = 1, and we have introduced new quaternionic quantities e'^, 63°. In 
all expressions a bar over a quaternion denotes its adjoint. 

All of the necessary functions of z entering the caloron formulas can be expressed in 
terms of these functions. From the equations that they satisfy (Eq. [2.11 , 2.10) we can 
easily deduce the general structure 

a a 

j;(^+D(+) + vI/-Z)(-)) (5.9) 



where the coefficients {Aa, etc) are quaternionic functions of space-time. 

For all the necessary calculations of adjoint modes one also needs the general form of 

duj: 

dcv = Y^i^TTiiz - Z'')^:Di+^ + + M/-5(-)) (5.10) 

a 

Finally, we need to know the integrals of the basis functions over z. We will need: 

It = / d.*,+t.„*j = ^ + n-) (5.") 

and 

t = J,..^r..K = (^) n°("->.) (5.12) 



- 12 - 



where we have introduced the quaternions 

p^" = i(^"±n,c7°n„) (5.13) 

The symbol ha stands for a hermitian unitary traceless matrix defined through the decom- 
position 

{x - X'')t = raha (5.14) 

where is the distance to the corresponding constituent monopole. The expressions also 
contain the function g: 

g{u) = uj sinh(M) (5.15) 

and its derivatives evaluated at the product of the mass and the distance. 

With the previous expressions one can compute the caloron vector potential as well as 
the adjoint modes, once the coefficients -Oi^^ and s')t'^ are known. These can be deduced 
from the equations that define u and lo (Eqs. p. 11 and 2.1C| ), but now the whole analytic 
structure reduces to a finite-dimensional linear problem in quaternions. Essentially, this 
follows from the matching at the edges of the intervals z = ^b\. For example, the absence 
of derivatives of delta functions in the equation for cj, implies that this function must me 
continuous at z = ±6i. If we note u;(5i) = Wi + iW2, then we can compute the coefficients 
of u and uj in terms of Wi by the continuity equations: 

W = e"^Z)i+) + e'^Di'^ (5.16) 
dW = -47r5ae"Di+) + e"'Si+^ + e'^S^-^ (5.17) 

The equation has been rewritten as a vector equation in terms of two (quaternionic) com- 
ponent column vectors W, e'"", ... . Notice that each equation is valid for both values of 
a. With ordinary vector space techniques one can solve for the coefficients: 

^(+) ^ m,ff(m,rj ^,^^ ^^ .^^^ 

rriara 



^(_) ^ -ihag{mara) 



mar, 



e^m (5.19) 



a' a 



^(+) ^ inag{mara) ^a^^^^, _ m,cosh(m,r,) ^,^^^^ ^^^^^^ 



maTa ' sinh^ (mara 

^(_) ^ -iha9{mara) _a^^^ ^ nia ^.^^ ^^ ^1) 

maVa sinh^ ( 

The coefficient Aa appearing in the expansion of u can be related to Dji~^ by the equation 
Mlo = u/F. Hence, we get 

Aa = 2iraFhaD\+^ (5.22) 
Combining all the previous formulas we arrive at 

uiajuja = W^CaaW + W^Ca&JW (5.23) 
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where 



_ 2Fmag{mara) ( . g'^imgrg) - I _ . g^(m^rj ^^^ , . 
^■^^ " - M ' 2m2r2 ^ + ' 2mara ' ^ ^ ^ 



vr 

and 



-iFg{maraj / , x_t + \ 
>Ca = eariQ -g[mara)el + (5.25) 

Up to this point all expressions seem to depend only on a single distance Tq. The 
mixing among the two coordinates and the relation between the constituent monopoles is 
hidden in the expression of W . The main equation satisfied by W is 

vr ^ ma 2 V -iT3 J 

From here one can solve for W . It is very easy to realise that W must be a linear combina- 
tion of the quaternion ir^j and unity. This follows from the equation oj = R^^q. Since q\s a, 
combination of these two quaternions and R commutes with quaternions (and is therefore 
real), this property extends to W . In summary, we have 

(y=fv(J^J (5.27) 

where V is a real 2x2 matrix, whose inverse is sum of contributions from the two constituent 
monopoles 

V-^ = P^ + Y,Ua (5.28) 

a 

with 

_ gimgrg) / cosh(mara) - cos(maXo) eaSin(maXo) \ ^g) 

2-Kma V eaSin(maa;o) cosh (mara )+ cos (ma^o) / 

An interesting relation between U and the vectors is given by 

vr rua 2tt \ —lUa / 

Now we have all the ingredients to calculate all the relevant quantities concerning 
calorons, including the adjoint zero-modes. For example, we can obtain the scalar function 
F 

For the vector potential, we can write the following expression: 

\ = ^W^a^diUi + U2)W+ h. c. (5.32) 

Finally, we proceed to the computation of the basic quantity which enters into the 
expression of the zero modes E^{x). This is obtained from Eq. 5.23| after multiplying by 
i and taking the hermitian part. The result is the sum of two terms. The first one has a 
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very transparent interpretation. To show this, one must first reahze that the traceless part 
(or hermitian part) of the quantity inside parenthesis in Eq. 5.24 is precisely E^^^{x — 



X°-;ma)/m'^, where E^^^{x — X°-;ma) is the gauge field of a BPS monopole of mass Ma 
centered at one of the constituent monopoles. This is sandwiched between the quaternion 
Qa = e°'^W and its adjoint. If we write 

Qa = e'^^W = \Qa\^i (5.33) 

where is a unitary matrix and we define 

\a = IQaP^ (5.34) 



vrma 

then we conclude that the first term in E'^{x) is given by 

\a^aEf'^{x-X''-ma)^i (5.35) 

which is just the field of a BPS monopole, gauge rotated and weighted by \a- The weight 
factors are positive and satisfy: 

Y.^a = {l-^) (5.36) 

a 

An explicit formula to compute them in terms of the matrix V is 



= 2 V'^^^-^^ + ^("^3 - X-,){Tt\V) - T,{V')) ] (5.37) 



The second piece contributing to Eq. 5.23| is also simplified if one realizes that 



Ca = -FdUa (5.38) 

Hence, one reaches to the following simple formula for the quantity i?^, representing the 
contribution of constituent monopole a to the field strength 

= ^g(^ara) ^tg^^BPS^^ _ x'^-ma) e^^W - —W^BUaoJw + h.c. (5.39) 
^ irma 2 

The time-like component (/x = 0) of the first term vanishes. For the second term we have 
£"0 = —Eq. This is easily concluded by realizing that 



dW = -Vd{Ui + U2)W (5.40) 
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